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. ' Abstract 

The asymptotics of the SU(2) 15j symbol are obtained using coherent states for the 
boundary data. The geometry of all non-suppressed boundary data is given. For some 
boundary data, the resulting formula is interpreted in terms of the Regge action of the 
geometry of a 4-simplex in 4-dimensional Euclidean space. This asymptotic formula can be 



Sh , used to derive and extend the asymptotics of the spin foam amplitudes for quantum gravity 



models. The relation of the SU(2) Ooguri model to these quantum gravity models and their 
continuum Lagrangians is discussed. 



£^ 1 Introduction 
O 

Q\ 

Spin foam models for quantum gravity are state sum models based on spin networks of various 
', kinds. In these models, the variables are representations and intertwiners of a group (or quantum 

group) labelling the edges and vertices of the spin networks. The amplitude for a simplex is 
defined by gluing these according to the combinatorics of the simplex, to make a closed spin 
network. The first model of this type was the Ponzano- Regge model, for three-dimensional space- 
times [PR68], where the 3-simplex amplitude is a 6j-symbol. Models in this spirit were developed 
for four-dimensional space-times in |BC00l IBTMl IKLPR081 IEPR081 lEPROTl IPer08l lLS08l IFK08] . 
These can be viewed as modifications or constrained versions of the much simpler Ooguri model 
|Oog92[ i ICY93j corresponding to a quantisation of BF theory. This topological model is based 
on the SU(2) 15j-symbol. 

A first insight into the geometric interpretation underlying these models can be seen in the 
large spin behaviour of the amplitude for a single 4-simplex. The intertwiner degrees of freedom 
can be identified with the shape space of tetrahedra with areas given by the spins [B98, BB99, 



IFKL091 ICF09] . R was shown in |BDF + 09a that if the data is such that these tetrahedra, as a 
geometric boundary, induce a non-degenerate Euclidean metric on the 4-simplex, the large spin 
asymptotics of the 4-simplex amplitude of the gravity models with Immirzi parameter [ELPR08, 
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IKPR081 IEPR071 IPer08l lLgQ8l IFK08] can be understood in terms of this geometry. In particular 
the phase part of the asymptotic behaviour is the Einstein action associated to this geometry 
in Regge calculus |Reg61j, as is expected in a discrete quantization of general relativity and as 
is the case, for example, in the Ponzano Regge model [PR68J. In this analysis, four terms were 
found to contribute to the asymptotics with phases proportional to the Einstein action. While 
it was expected to find two contributions from the two orientations associated to this boundary 
data, the existence of additional terms came as a surprise. 

In this paper we apply and extend the methods of [BD F + 09a to study the asymptotics of the 
SU(2) 15j-symbol underlying the new models. We find that the asymptotics are governed by the 
expected classical solutions of the topological theory. We find however that for some boundary 
data these solutions retain an interpretation in terms of 4-dimensional Euclidean geometry and 
the phase part of the asymptotics of the 15j-symbol is given by the Regge action. This provides 
a proof for a part of the asymptotic formula for the the Lorentzian model BDF + 09b] that was 
missing. The obtained asymptotics for the 15j-symbol are then used to derive the asymptotic 
formulae for the four-simplex amplitudes of the new Euclidean gravity models expressed as 
rescaled squares of the 15j-symbol. This provides a nice explanation for the appearance of 



the weird terms in BDF + 09a . A further explanation for the non-geometrical terms in the 
asymptotics of the gravity models is explained by the fact that the gravity models appear to 
contain a single copy of the Ooguri model as a subset of the configurations. This interference 
between the topological Ooguri model and the gravity configurations has some implications for 
the continuum limit of the gravity state sum models, as discussed in section HJ Finally, there 
is a heuristic interpretation of our results in terms of the corresponding Lagrangian continuum 
field theories. 



2 SU(2) 15j-symbol asymptotics 

The 15j-symbol is determined by data associated to the boundary of a 4-simplex. This data is a 
spin k in N/2 for each triangle and a four-valent intertwining operator for each tetrahedron. By 
introducing coherent states for the spins k, one can represent each intertwiner in the coherent 
state basis (LS07J by SU(2)-averaging the tensor product of four coherent states. The notation 
used is as follows. A coherent state for the spin ^ representation is labelled by a unit vector n 
in M 3 , and is denoted |n); the phase of this coherent state is arbitrary at present and will be 
specified later. The tetrahedra in the 4-simplex are labelled with a = 1, 2, . . . , 5. The spins label 
the triangles of the 4-simplex and so are specified by k ab = k ba , for each a ^ b. 
The intertwiner for the a-th tetrahedron is given by the formula 

4>a= [ (g)X|n a6 > 2 ^, 

Jxesv(2) h+a 

and the boundary state for the whole simplex is 

V = (g)Vv 

a 

The data {k ab , n ab } specifying this boundary state up to phase is called the boundary data. The 
phase choice will be discussed later. 

The 15 j symbol in the coherent state basis is given by the following integral formula [LS07], 
|BDF+09a| 

15i(fc,n) = (-l)* / TT dX a TT (Jn ab \XlX b \n ba ) 2k « b . (1) 

M2)5 atts l<a<b<5 
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Here, f is the Hermitian conjugate, J is the standard anti-linear SU(2) operator and the notation 
(Jn a t| stands for the covector dual to J \n a b)- Finally, the sign factor ( — l) x is determined by the 
graphical calculus relating the 15 j spin network diagram to the above evaluation. The formula 
is a linear function of the boundary state ip and is determined by contractions, according to the 
combinatorics of the diagram, using the bilinear inner product on spin k representations. 



2.1 Asymptotic problem and critical point equations 

The above formula is an integral expression in exponential form and so the asymptotic limit, 
where all spins are simultaneously rescaled, k ab — > Xk ab , and taken to be large (A — > oo), can be 
analysed with stationary phase methods. 

The action for the asymptotic problem yields 

%n) [X] = 2k ab In ( Jn a6 | X\X h |n 6a ) . (2) 

a<b 

It enjoys a global SU(2) continuous symmetry and a discrete ± symmetry at each vertex a, given 
by the formula 

X' a = e a YX a , (3) 

with Y £ SU(2) and e a = ±1. 

The critical points of this action satisfy the same equations as in [BD F+09a| . The stationarity 
with respect to the group variables implies closure for each tetrahedron a 

^2 KbTO-ab = 0, (4) 

b:b^a 

and the critical points are those which also satisfy 

X b n ba = -X a ii ab (5) 

for each triangle ab, a / b. 



2.2 Solutions to the critical point equations 

In this section, we firstly show that the solutions to the above critical point equations admit 
an interpretation in terms of BF theory. Then, we study the classification problem of these 
solutions for various classes of boundary data. 



2.2.1 Geometry of the solutions: relation to BF theory 

The solutions to the critical point equations @ and ([5]) can be interpreted as the solutions of 
a four-dimensional BF theory with group SU(2) discretised on a 4-simplex. The fundamental 
fields of the theory are an SU(2) connection A and an su(2)-valued two-form B. The equations 
of motion of the theory are solved by covariantly closed B fields and flat connections: 

d A B = 0, and F A = 0, 

where Fa is the curvature of A. Note that since flat connections are locally pure gauge, the first 
equations reduces locally to dB = 0. 

Evaluating the BF action requires an orientation of the 4-manifold. In the discrete setting, 
an orientation of a 4-simplex establishes a coherent rule for associating an orientation to the 
triangles. Triangle t ab is defined to be the intersection of tetrahedra a and b with the orientation 
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induced from the boundary of tetrahedron a, which itself inherits an orientation from the 4- 
simplex. Thus, we have that i& a = — t a b- 

The analogy between BF theory and the asymptotic data is apparent by defining the fol- 
lowing variables 

b a & = kabXaYlab. 

The solutions to the critical point equations can then be parameterised by the X a and h a b- 
These are the discrete connection and 5-field variables, respectively. This interpretation follows 
from the critical point equations (£[]) and ([5]), which expressed in terms of the h a b, read 

b a6 = 0, b ab = -h ba - (6) 

b:b^a 

A such set of twenty three-dimensional vectors h a b satisfying the above two vector equations 
determines a geometrical structure called a vector geometry. Accordingly, a vector geometry is 
parametrised by 20 x 3 — (4 x 3 + 10 x 3) =18 numbers. This data determines a constant 5-field 
on the 4-simplex, namely the 2-form valued in the Lie algebra su(2) with coefficients which are 
constant in the standard linear coordinates on the simplex, satisfying 

f B = h ab . 

•> tab 

The closure equation Yl b a 6 = is now interpreted as Stokes' theorem for the B field around 
the boundary of a tetrahedron, that is, the first equation of motion of BF theory dB = 0. 

The proof that the existence of the constant S-field is equivalent to ([6|) is as follows. Pick a 
distinguished vertex of the 4-simplex and denote the four edge vectors pointing away from that 
vertex by £i, ... ,£4. A basis for the space of bivectors A 2 (M 4 ) is given by £ p A £ ? . Therefore the 
components of B with respect to this basis are -B(£ p A £ 9 ) which are just the integrals of B on 
the triangles meeting the distinguished vertex. Using the closure relation shows that these are 
all the independent variables in a solution to ([6]). 

The X a variables are a discrete version of the connection. This is because the gluing of 
two 4-simplexes proceeds by the identification of the n a b variables on a common tetrahedron. In 
terms of the h a b variables, this means the gluing takes place after the action of the corresponding 
X a for the tetrahedron. The X a is therefore a parallel transport operator for half the dual edge 
between the two 4-simplexes. As the boundary data has not been varied in the critical equations, 
the flatness equation of motion for this connection is not yet enforced because it is invisible at 
the level of a single 4-simplex. This equation appears only when 4-simplexes are glued all the 
way around a triangle. 

2.2.2 Classification of the solutions 

The existence and a partial classification of the solutions to the critical point equations (JH) and 
© was established in |BDF + 09a] and depends on the boundary data. Two solutions X a , X' a 
are said to be equivalent if they are related by the symmetry operation ([3|) and inequivalent 
otherwise. In theorem [1] the results of |BDF + 09a] are extended by studying the existence and 
classification of these equivalence classes of solutions as a function of the boundary data. 

General boundaries. For general boundary data, the solutions to the vector geometry equa- 
tions ([B]), for fixed |b a b|, were analyzed by Barrett and Steele |BS02j. The vectors for a given 
tetrahedron describe a quadrilateral in M 3 . Thus there are no such solutions unless the |b &| 
satisfy the inequalities of a quadrilateral (i.e., each one can be no larger than the sum of the 
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other three). Given a solution to ©, one can glue up an open cubical box starting at any 
quadrilateral. The equations imply that the 5 quadrilaterals glue up as 5 sides of a box and the 
open side of the box is a parallelogram. Conversely any such box geometry determines solutions 
to ([5]). Barrett and Steele showed that for a given set of k a b, determining the edge lengths of 
the box, there is a 5-dimensional set of such box geometries. 



Non-degenerate boundaries. There is an important subset of boundaries where the above 
results also apply. Boundary data such that for each tetrahedron a, the face vectors n a b, for 
fixed a and varying b, span a three-dimensional space is called a non- degenerate boundary data. 
In this case, if the four vectors n a b satisfy the closure condition (jl]) they specify an embedding 
of the tetrahedron in three-dimensional Euclidean space, such that the vectors are the outward 
face normals and the k a b are the areas. In this way, each tetrahedron inherits a metric and an 
orientation but the metrics and orientations of different tetrahedra do not necessarily match. 
For non-matching boundary data, the classification established in BDF + 09a] states that there 
is, up to equivalence, one or no solution to the critical point equations. 



Geometric boundaries. Non-degenerate boundary data for the whole 4-simplex is said to 
be geometric or Regge-like if the closure constraint (J3]) is satisfied and if the individual tetra- 
hedron metrics and orientations glue together consistently to form an oriented Regge-calculus 
3-geometry for the boundary of the 4-simplex. This is the requirement that the induced metrics 
on the triangles agree for both of the tetrahedra sharing any given triangle, and the induced 
orientations are opposite. This boundary data satisfies the gluing constraints discussed by Dit- 
trich, Ryan and Speziale in [DR08, DS08J, which were found there to be necessary to ensure 
geometricity. 

For geometric boundary data, there exists a unique (up to a Z2 lift ambiguity discussed in 
BDF + 09aj ) set of ten SU(2) elements g a b = which glue together the oriented geometric 
tetrahedra of the boundary, defined as follows. The embedding of tetrahedron a in Euclidean 
space determines a position vector for vertex v. Suppose v and w are vertices which appear 
in both tetrahedron a and tetrahedron b. Then gb a is a spin lift of the unique rotation which 
maps — x^ to x^ — x™, for each possible choice of v and w, and maps the outward normal to 
one tetrahedron to the inward normal to the other, 

9ba n a b = -n 6a . (7) 

This gluing map can also be interpreted as the parallel transport operator gb a '■ T a — > T& mapping 
the tangent space T a of the a-th tetrahedron to the tangent space Tb of the b-th. tetrahedron_| 
and describing the change of reference frame from a to b. 

From this data, one can make a canonical choice of phase for the boundary state ip by picking 
the phases of the coherent states such that 

|n& a ) = 9baJ\^-ab)- (8) 

The boundary state ip with this choice of phase is called a Regge state. 

Geometric boundaries of a 4-simplex fall in three geometrical classes, depending on whether 
it is the data of a 4-simplex with 4d Lorentzian (with space-like tetrahedra), 4d Euclidean or 
3d Euclidean geometry. This follows from the fact that a metric for the 4-simplex is determined 

1 Note that to account for this our convention here differs by g a t — * gt a from [BDF + 09a] , The vector n a (, is 
interpreted to be in tetrahedron a. With this convention g^a goes from a to b which is notationally more convenient 
and in accordance with the parallel transport conventions. 
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uniquely by the edge lengths. Since the three-dimensional hypersurfaces containing the tetrahe- 
dra all have signature + + +, the signature for this 4- metric can only be — h ++, + + ++ or 
+ ++. In the last case, the metric is a pull-back of the Euclidean metric on M 3 by a projection. 
To complete the description of geometric boundary data initiated in [BDF + 09a] and [BDF + 09b] , 
where the boundary data for Lorentzian and Euclidean 4-simplexes were discussed, we briefly 
comment on the boundary data of a 4-simplex in M 3 . 

This case is a 4-simplex with a metric of signature H — I — |— and has the same metric geometry 
as a linear immersion of the simplex into M 3 . Thus it follows that the map from the 4-simplex 
boundary to M 3 has to be orientation-reversing on some simplexes. A geometric way to see this 
is to imagine a 4-simplex in 4d and gradually flatten it. Some of the outward normals to the 
tetrahedral faces of the 4-simplex will point into the northern hemisphere, while others go to 
the southern hemisphere. As the boundary orientation plus the outward-facing normal agree 
with the standard orientation on M 4 , the tetrahedra with southward-pointing normals will be 
orientation-reversed with respect to ones with northward-pointing normals. Let us suppose that 
the standard orientation in M 3 agrees with the tetrahedra with southward-pointing normals. 
Then applying the inversion map M 3 — > M 3 

/:xh —x, 

to all the tetrahedra with northward-pointing normals gives all of the boundary tetrahedra 
embedded in 3d Euclidean space with orientation-preserving maps. This embedding of the 
tetrahedra in M 3 defines a set of outward normal vectors to the triangular faces of the tetrahedra 
n a fc, and hence a set of boundary data for this geometry. Moreover the gluing maps gb a have 
a straightforward description with this representation. If two tetrahedra are both northern 
hemisphere or both southern hemisphere, then <#, a is a spin lift of the identity rotation. If one 
is northern and one southern, then g^ a is a spin lift of the rotation by n about the normal to 
the common triangle in each tetrahedron, n a b = — n& a . A solution to the critical equations © 
is given by X a = 1. The boundary data just described is not the only possible set for this 
geometry. All other sets of boundary data for this geometry are obtained by independently 
rotating each tetrahedron. This defines boundary data in the 3d Euclidean case. 

We are now ready to extend the classification and geometry theorems of [BDF + 09a]. 

Theorem 1. Let the boundary data of a 4-simplex be geometric. Then the critical point equa- 
tions (J5J) have 2 equivalence classes of solutions in the case of 4d Euclidean boundary data, 1 
equivalence class in the case of 3d Euclidean boundary data, and no solutions in the case of 4d 
Lorentzian boundary data. 

Proof. Part of this theorem is proved in [BDF + 09aj . Specifically, it is shown there that there 
are 2 equivalence classes of solutions if any only if the boundary data is 4d Euclidean. 

The remainder of the theorem is proved in the rest of this section by giving a general proce- 
dure that operates on a given solution {X a } to give a second solution {X' a }. The construction 
proceeds as follows. 

A solution {Xa} to the critical point equations can be written alternatively using the variables 

X a b = X a 1 Xb, 

which satisfy 

X a bXb c X ca = 1, 

for all a, b, c. The critical point equations ([5|) then become 

X a b n ba = ~ n ab- 
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The involution which generates a second solution {X ( '-^a^b °f twenty group elements is con- 
structed from the gluing maps and the first solution as follows 

X'ab = 9ab Xba 9ab- (9) 

In fact, since X' ba = (X'^) -1 by construction, we have ten group elements {X' ab ] a< b satisfying 
the equations 

X'at, n& a = — Il a b. 

The set of group elements {X' b } a< i) are however not solutions to the critical point equations 
unless each X' ab can be expressed as X' ab = (X'^^X^ for some X' a and X' b in SU(2). This is 
equivalent to proving a cocycle condition for each face abc of the complex dual to the 4-simplex, 

X 'ab X bc X 'ca = ^ 

because the fundamental group of the three-sphere is trivial, and so all flat connections are pure 
gauge. To show that the above cocycle condition is indeed satisfied we need the following lemma. 

Lemma 1. Let A,B,C be three SU(2) matrices corresponding to three SO(3) elements with 
rotation axes lying in the same plane P. Let R be an SU(2) element associated to a rotation in 
the plane P. If ABC = R then CPA = Br 1 . 

Proof. Let X be a 7r-rotation in the plane P. Conjugating both sides of the equality 
ABC = R by X leads to 

XAX' 1 XBX^XCX- 1 = XRX- 1 . 
Now, it is immediate to see that XDX' 1 = D' 1 , for D = A, B, C, and XRX' 1 = R. □ 
The above lemma is applied to our framework by defining the following three SU(2) matrices 

A = g a bXba, B = g a bgb c X c b 9ba, and C = 9ab9bc9caX a c 9cb9ba- (10) 

These matrices satisfy CBA = D a b c , with B> a bc '■= 9ab9bc9ca- Furthermore, using ([5]) and (J7J), it 
is straightforward to check that the following three vectors 

s A = n ab , s B = g a bn bc , and s c = g a b9bcn C a, 

are stabilised by A, B and C respectively. This implies that the axes of the rotations associated 
to A, B and C are given by s^, s# and sc respectively. 

Now, recall that the gb a map corresponds to the parallel transport operator from tetrahedron 
a to tetrahedron b ; it encodes the change of reference frame from a to b. This leads to two 
conclusions. Firstly, the SU(2) element P> abc = g a b9bc9ca, mapping the tangent space of the 
tetrahedron a to itself, corresponds to a rotation in the plane P orthogonal to the edge (abc) 
(hinge) shared by the three triangles ab, be and ca. This rotation stabilises the corresponding 
edge vector e a = — x^ 1 , constructed from the position vectors of the vertices v,w bounding 
(abc) in the frame associated to the simplex a, and the associated angle is the deficit angle 
around this edge. Secondly, because the normals n a b, rift c and n ca correspond to three triangles 
sharing the abc edge, they can all be defined in terms of the corresponding edge vector. Indeed, 
the normals can be constructed as follows 

n a b = e a x u a6 , n bc = e b x u bc , and n ca = e c x u ca . 
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where, the symbol 'x' denotes the 3d cross product, and u ab , u bc , u cc i are a second, distinct 
edge vector in each one of the triangles ab, be and ca respectively, which is chosen such that 
n ab is pointing towards the exterior of simplex a. Now, e a , e b and e c are related by the parallel 
transport operator and we have that e b = gba^a, arid e c = g cb g ba e a - Therefore, one can readily 
see that 

n a fo • e a = 0, n bc ■ (g ba e a ) = 0, and n ca ■ (g cb g ba e a ) = 0, 

where the dot '•' is the 3d Euclidean inner product. Using the invariance of the associated 
metric, that is, (gx) • y = x • (<7 _1 y)) we can therefore see that s^, s# , and sq, the axes of 
rotation of the three matrices A, B and C, all lie in the plane . 

We can therefore apply lemma [T] to the present context with the role of the matrix R played 
by D abc and the plane P = is given by the orthogonal complement of the edge (abc). This 
leads to the conclusion that ABC = which reads, for each edge (abc), 

X 'ab X bc X 'ca = I- (11) 

Therefore, for each triangle (ab), there exits two SU(2) element X' a and X' b such that 

X 'ab = ( X 'a) X h 

and we have generated a second solution {X' a } to the critical point equations. There are now 
two cases to consider. The two solutions X and X' are either inequivalent or equivalent. 

• If the two solutions are inequivalent, the reconstruction theorem of [BDF + 09a] which states 
that the vector geometries determined by {X a } and {X' a } recombine into a single bivector 
geometry applies and the boundary data is 4d Euclidean. 

• It is also possible that the two solutions X and X' are in fact equivalent, that is, they are 
related by the symmetries ©, which implies that X' ab = e a e b X ab . Now if the boundary 
data for each tetrahedron is rotated so that the original solution is X a = 1 for all a, then 
([9]) implies that the gluing map satisfies 

alb = ±1 

for each a ^ b. This means that the corresponding rotation squares to the identity in 
SO(3) and is thus either the identity element or a rotation by n about the axis n ab . To 
discuss the pattern of these two possibilities, define a parameter a ab to take the value — 1 
for a rotation by ir and +1 for the identity rotation. 

To show that this description of the boundary data and gluing maps g ba coincides with 
that of a 3d Euclidean geometry, it is necessary to show that the tetrahedra fall into two 
classes, with a ab = — 1 for a transition between the two classes and a ab = +1 if both 
tetrahedra lie in the same class. Consider three distinct tetrahedra, labelled with a, b and 
c. These tetrahedra have a common edge, which is represented by the displacement vector 
e a in tetrahedron a (v and w labelling its vertices, as above). Then 

gaegebgba^a — 6 a . 

But g ba e a = a ab e b , etc., so that 

O'abO'bcO'ca = 1- 

It follows that 

&ab = 0- a a b , 

with a a = dbl. Accordingly, the tetrahedra fall into two groups, those with a a = +1 and 
those with a a = —1, and g ba is a rotation by tt only when gluing one of one group with 
one of the other. This shows that the geometry is exactly that of the 3d Euclidean case. 
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To complete the proof of theorem [U the argument is as follows. The explicit description of 
the 3d Euclidean boundary data shows that at least one equivalence class of solutions of ([5]) 
exists for this case. There cannot be more than one because the part of the theorem proved 
in |BDF + 09a would then show that boundary data is that of a non-degenerate 4d Euclidean 



4-simplex. 

Conversely, if a solution to ([5]) exists, then we obtain another one by applying the involution. 
By the explicit analysis of the case where the involution yields a symmetry-related solution 
we saw that this case corresponds to the boundary data of a 3d Euclidean 4-simplex. If it is 
inequivalent, the part of the theorem proved in [BDF + 09a applies again, and it is 4d Euclidean. 
It follows immediately that for the Lorentzian case there are no solutions to these equations. □ 

2.3 The action at the critical points 

In this section, we evaluate the action on solutions to the critical point equations and derive 
asymptotic formulae for the 15j symbol in the case of non-degenerate boundary data. 

2.3.1 Non-geometric boundary data 

Lifting the critical point equation to coherent state space produces a phase 

(X a )iX b \n ba ) = e** J\n ab ), (12) 

for all a ^ b, from which follows that <p ab = <p ba . This phase is the evaluation of the action as of 
course 

(Jn ab \(X a )^X b \n ba } 2k ^ =e 2ik ^\ 
and thus the action ([2]) evaluates to 

S(k,n)\crit = 2i k ab <pab- (13) 
a<b 

When scaling up the k this will provide exactly the phase of the oscillating part of the asymptotic 
formula. Note that the phase <j) ab depends on the arbitrary phase of the coherent states. While 
this arbitrary phase cancels in the full state sum we need to make a choice when decomposing 
the state sum into individual weights. 

Under the assumption that the boundary data is non-degenerate (although non-geometric), 
we have seen that there is at most one solution to the critical point equations (JSJ up to equiv- 
alence. Thus, we can change the \n ab ) to e l ^ ab \ rv ab ) for a < b, and set the action at the critical 
point to zero. This choice of course depends on the full set {k ab , n ab }. Using this phase choice 
we can now easily evaluate the asymptotic behaviour of the 15j symbol 

l^W-l)*^)'^^, (14 ) 
where H is the Hessian of the action ([2]) evaluated on the critical points and the square root is 



defined as in BDF+09a . 



2.3.2 Geometric boundary data 

If the boundary data defines a positive definite metric on the boundary of the 4-simplex we have 
a geometrically distinguished boundary state, the Regge state, defined by ([5]). We can then look 
at the three cases of the classification theorem in turn with the boundary state given by a Regge 
state. 
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Two inequivalent solutions. Given two inequivalent solutions b^ b = k a bX^n a b, we know 
from the reconstruction theorem in p3 DF + 09a] that we have a pair of inversion-related non 



degenerate 4-simplices a with bivectors given by their self-dual and anti-self-dual parts by 
Bah(?) = (&o&> b~ b ), together with a pair with parity-related geometry with bivectors given by 
B a b{P&) = — (b~ b ,b^ b ). This notation fixes our choice of ± labels corresponding to the self-dual 
and anti-self-dual components of the bivectors for the two solutions. 

From the lift (|12p we now obtain two set of angles (fy^ h . By substituting © and ([8]) into (|12p 
we have immediately that 

„id>\ „—i<t>~u 

The more detailed analysis in |BDF + 09a] shows that these angles are, up to factors of tt that 
cancel once exponentiated, equal to plus or minus half the dihedral angle Q a b of the triangle 
shared by tetrahedra a and b 

1 

-( 

"2 

As the k a b give the areas of the geometry the action thus evaluates to the Regge action Sr(ct) 
for the 4-simplex a determined by the boundary data 

%,n)l± = ±*X^ kab ® ab = ±iS R (a). ( 15 ) 

a<b 



4 = 



The full asymptotic behaviour is then given by: 

/2vr\ 6 2 4 / e iXS *W e -i\S R (a)\ 

m », n , ~ (_!)« (_) _ + 7sr _j , (16) 

with H± the Hessian of the action evaluated at b± respectively. 

Note that as we have two terms now an arbitrary phase change in the coherent states could 
at most set one of these phases to zero, (ft + — <p~ = O is independent of the arbitrary phase. The 
geometric phase choice is significantly stronger than that in the non-geometric sector as it can 
be applied directly to the boundary data, without having to solve the critical point equations. 



Single solution. In this case we have seen that we have the geometry of a degenerate 4-simplex 
in M 3 and a solution is given by the identity rotations. Then the phases (j> a b evaluate to the phase 
of gba- The action is a degenerate case of (|15p. with @ a b = or tt. As per the discussion of 
the flat geometry we have Q a b = tt exactly if a and b are oppositely oriented. Therefore in the 
action tt is multiplied by a sum of spins belonging to a closed surface and these always sum to 
an integer. Since the action is a multiple of tt, the ± sign is irrelevant, and e lSR ^ = ±1. The 
full asymptotic behaviour is given by this sign times (|14p . 



No solutions. This situation occurs when the boundary geometry is that of a 4-simplex with 
Lorentzian metric. In this case, the 15j symbol goes to zero asymptotically faster than any 
polynomial of A. 



3 Quantum gravity asymptotics from the 15j symbol 

In this section, we show how the previous results relate to models of Euclidean quantum gravity. 
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3.1 Quantum gravity models as constrained squares of the Ooguri model 

All the new quantum gravity models are closely related to the Ooguri state sum which is con- 
structed using the SU(2) 15j symbols. In fact most of them can be written as constrained 
versions of the square of the Ooguri model in the coherent state basis. The state sum model is 
given by identifying the boundary data on each 4-simplex appropriately and integrating over it. 
The boundary data for the gravity model on a 4-simplex is a set of unit vectors n ab and a set of 
spins j a b, in other words exactly the same data as for an SU(2) 15j symbol. The precise details 
of the gluing of simplexes are beyond the scope of this paper, but some further remarks about 
the full state sum are made in section 21 

There are two types of model given by products of 15j symbols. The first has 4-simplex 
amplitude given by rescaled squares of the 15j symbol 

n) = 15j(j ab , n ab ) x 15j(c 7 j ab , n ab ). (17) 

using the parameter c 7 = . 

This gives both the EPRL and FK models for parameter 7 < 1, and on setting 7 = (i.e., 
c 7 = 1), the EPR model. 

The second type of model has 4-simplex amplitude given by 

Z cc (j,n) = 15j(j ab ,n ab ) x 15j(c 7 j a ft, n a6 ). (18) 

and includes the FK model with parameter 7 > 1, and for 7 = 0, the FK model without a 
parameter. The only model that can not be written in this form is the EPRL model with 7 > 1 
which takes the form Z cc asymptotically. 

Note that these models are usually written as function of k such that j = ^-ip-k and f = 

^Ip' k. The possible values of 7 and j, j' and thus k are restricted as the state sum is identically 
zero unless c 7 j takes on half integer values. 

This formulation was the way the FK model [FK08J was initially written down, and for the 
EPRL model this formulation was developed in [LS071 lBDF+09^] . 

3.2 Quantum gravity 4-simplex asymptotics from the SU(2) 15j symbol 
The 4-simplex amplitudes for the models without complex conjugation are then given by 

Z(j, n) = (-1)*' f TT dX a dX' a ( Jn ab \ X\X b |n 6a ) 2 ^ (Jn ab \ X ,] a X' b |n 6a ) 2c ^ . (19) 

For the models with complex conjugation, the 4-simplex amplitudes read 

Z cc (j,n) = (-I)*' f TT dX a dX' a (Jn ab \ X^X, \n ba ) 2 ^" (Jn ab \ X'\X' h |n fea ) ^ . (20) 

^SU(2)io a ^ 5 

The corresponding actions can be written easily in terms of the 15 j action ([2]): 

% n) [X,X'] = 5 - n) [X] + c 7 % n) [X'], 

S^[X,X'] = S Utn) [X] + (h S Utn) [X']. (21) 

The equations governing the stationary phase analysis of all these models are identical and are 
simply (UD and two copies of (0): 

jab^ab = 0, 

X b n ba = —X a n ab , 

X' b n ba = -X' a ii ab . (22) 



11 



We can thus directly apply the geometric discussion above to the solutions to these equations in 
terms of b ab = j a bX a n ab . In particular it is clear that given any set {b a b,X a } solving the SU(2) 
equations we obtain a solution of the equations for the new models simply by setting X' a = X a . 
Generally, the asymptotic behaviour can be obtained by simply multiplying the asymptotics of 
the two SU(2) factors. 

3.2.1 Non-geometric boundary data 

We obtain two angles, 4> ab and <j>' b from the lift of (|22p . These combine to give the summands 
in the action of the models without complex conjugation 

(Jn afe | (X a )^X b \n ba ) 2 ^ (Jn ab \ (X'J X' b \n ba ) 2c -<^ = e i2 ^> +c -<^\ (23) 

Analogously for the models with conjugation the phase is (p ab — c^4> ab - Thus the actions evaluate 
to 

S Jj,n) l c «i = 2i ^2 jab(4>ab + C^(j)' ab ), 
a<b 

SJj%)\crit = 2i^2jab(4>ab ~ C~,(j)' ab ). (24) 
a<b 

Restricting to non-geometric boundary data which are non-degenerate, we know that there 
is one or no solution to the critical point equations ([5]), up to equivalence. Thus we have that 
{Xa} ~ {X' a } and therefore that 4> ab = (j)' ab . Thus, just as in the case of a single SU(2) we can 
now set the phase in \n ab ) such that it cancels the (j) ab to obtain the final asymptotic formula by 
multiplying two copies of (Q3I scaled by c 7 

ZW.nW-l^f)' 2 ^;^^, (25) 

where H is the same Hessian of the action ([2]) evaluated on the critical points with spin given 
by j. As the Hessian is linear in the j and 12 dimensional, the determinant evaluated at the 
critical point f = c^j differs by c 7 12 from that evaluated at j. This accounts for the factor c 7 6 . 

3.2.2 Geometric boundary data 

If the boundary data defines a positive definite 3-metric on the boundary of the 4-simplex, we 
again chose the boundary state to be a Regge state ©, and the classification proceeds in the 
same way as before. 

Two inequivalent solutions. As before given two inequivalent solutions b b = j ab Xaii ab they 
necessarily correspond to the self-dual and anti self-dual parts of geometric 4-simplex bivectors 
given by B ab (a) = (6+ 6 ,6^ 6 ) and B ab (Pa) = -(b~ b ,b+ b ). 

These combine to four solutions to the stationary and critical point equations (|22p : 

(6, a," 1 b') G {(&+ b + ), (6+ 6"), (&-,&+), (6", 6-)}. 

The ++ and solutions are analogous to the solutions for non- geometric boundary data as 

the solution is just the double of a pure SU(2) solution. With the Regge boundary state the 
phase in (fl~2|) and thus the phases (j), (j)' in ([23]) evaluate to 4> a b\b ± = =t^©ab- 
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The full asymptotics is then again given in terms of the Regge action Sr(ct) for the 4-simplex 
a determined by the boundary data by a sum over the signs of the 4 critical points {b^^b^ 1 ). 
That is, with e, e' = ±1, 

/2vr\ 12 2 8 ^ i\{e+e'c^S R (a) 
Z«,n)~(-1)*( T ) _ i £- 7 _ g __, ( 26) 

for the models without complex conjugation and by 

Z~(A;,n) -(-!)* (fl) * £ 1 (27) 

V A ; (47r) lb ^e^det^det^/ 

for those with. As the amplitudes leading to these are rescaled squares of the SU(2) 15j symbol 
this is of course just the rescaled square of the SU(2) asymptotics (fl~6j) . 

We can rewrite the actions further by expressing them in terms of k a b = r^x jab- Writing 
S tt / for 5*| (5 e ,6 £ ) we obtain 



S] t/ = ( e (l+ 7 )+ e '|l- 71)^^6^, 

a<b 

Sl? c = ( e (l + 7 )- e '|l- 7|)jE^ ^ 



2 



and thus 



gr>x = ggn* = ((£ + e , h + {e _ e / )} 1 ^ 



2 

a<6 



57 >l,cc = 57 <1 = ((£ _ £ , )7 + (£ + e/)) 1 £ ^ 0afe) (2g) 

a<6 



Hence, we have in particular 



SWlllee' — 5'^ < |ee' 7=0 — 2 ^ ] ^afeQqbi 



and 

e — e' 

SWIee' = S |ee' 7=0 = 7^ ^ ] kqb®ab- 

This is the form in which the actions appeared in the literature so far, where the state sums 
were expressed in terms of k. 

Single solution. As in the SU(2) BF case the action is zero or ir and the amplitude evaluates 
to (|25|) times a sign. 



No solutions. In this case the boundary geometry is that of a 4-simplex with Lorentzian 
metric and we again have that the amplitude is suppressed exponentially for large spins. 
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4 Discussion 



The SU(2) 15j-symbol is the building block for the Ooguri model which is a topological model 
based on flat SU(2) connections. Although it is much simpler than a quantum theory of grav- 
ity, the results of this paper have demonstrated that understanding the 15j-symbol and its 
asymptotics is the key to understanding the more complicated gravity models. Curiously, the 
15j-symbol exhibits some gravity-like features itself. 

There is a heuristic dictionary which translates between the discrete variables and discrete 
action entering into asymptotics formulae, and the traditional continuum fields and their actions 
constructed from Lagrangians. This dictionary has been started in the main text, where the 
interpretation of the general asymptotic data for an SU(2) 15j-symbol is interpreted as discrete 
B and A fields. This analogy goes further as it is possible to construct the Ooguri model from a 
BF Lagrangian, in much the same way as the Ponzano-Regge model is related to a BF theory 
in three dimensions. In fact, a configuration of the asymptotic data for an SU(2) 15j'-symbol is a 
vector geometry which gives a constant su(2)-valued B field on the 4-simplex and a discrete SU(2) 
connection determined by the X matrices. In the context of a state sum model on a manifold 
it is the sum over the spin labels k on each triangle that enforces the flatness of the holonomy 



around the triangle Oog92 . This flatness is not visible by examining a single 4-simplex. 

The analogy goes further. A B field is 'geometric' if it is the self-dual part of e A e for some 
non-degenerate frame field e 1 

B l = Pjje 1 A e J , (29) 

with Pjj projecting antisymmetric 4x4 matrices to their self-dual part, indexed by i. If B here 
is the constant su(2)-valued 2-form on a 4-simplex determined by a critical point, it satisfies this 
equation precisely when the boundary data is geometric. This is because a geometric boundary 
induces a constant metric on the 4-simplex, for which a frame field exists. One can see that 
there are five constraints involved in these equations. The constant B field has 18 components, 
whereas the geometric ones form a space of dimension 13. (The e field has 16 components but 
three-parameter's worth are identified in (|29p . due to the invariance of the self-dual part under 
anti-self-dual rotations.) 

A further confirmation of this interpretation is given by the continuum analogue of ([9]). In 
this equation, X represents the self-dual part of the connection uj + , X' the anti-self-dual part 
cj_ and g the Levi-Civita connection on the boundary ujq. Linearising (J9j) leads to 

u>- = 2u)q — u+, 

which is the right relation for these connection forms restricted to the boundary if the two-form 
B is geometric. 

This dictionary also gives a heuristic explanation for the appearance of the Regge action for 
geometric boundary data in the 15j-symbol. This is because in the continuum picture, the SU(2) 
BF action is exactly the action for self-dual gravity on substituting with ([29]) . as discovered by 
Plebanski. 

Hence the restriction to geometric or Regge-like boundary data is the precise counterpart 
in the discrete setting of equation (|29[) in a continuum Lagrangian. This suggests that it is 
possible to define a discrete version of Plebanski's gravity, in which constraints are added to a 
BF Lagrangian to enforce (|29p . The quantisation with a state sum model is given by the usual 
formula (here somewhat simplified) 



Z PL = J ^2 II 15j-symbol 

{ n } k 4-simplexes 
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but with the integral over the coherent states restricted to ones that specify geometric bound- 
ary data on each 4-simplex. Whether such a model is simple enough or useful remains to be 
seen. More generally, one could seek modifications of the Ooguri model so that the geometric 
configurations dominate and the Plebanski Lagrangian, and hence general relativity, is evident 
in an asymptotic regime. Such a modification would of course have to break the reduction of 
the Ooguri model to flat connections. 

The apparent paradox that the topologically trivial Ooguri model contains gravity in its 
asymptotics is resolved by the fact the sum over the spins in the Ooguri state sum model will 
always enforce flatness. Thus it is possible to have geometric boundary data for an arbitrary 
triangulation of a 4-ball but the asymptotics will always be determined by the Einstein action 
for flat metrics, as in the three-dimensional case [DGH09J. 

Moving now to the discussion of the gravity state sum models based on SU(2) x SU(2), the 
results here demonstrate that the asymptotics of these models can be explained by a suitable 
squaring of the SU(2) 15j symbol asymptotics. This gives a very nice explanation of the four 
terms in the asymptotics found for the EPRL/FK models. These are just obtained from squaring 
an asymptotic formula with two terms. In particular this demonstrates that the two of these 
terms that were unexpected (called the 'weird terms' in BDF + 09a] ) are definitely required. 



However the analysis of this result, which was the original motivation for the paper, delivered 
also another insight into the gravity models: the asymptotic formula always contains the terms 
which correspond to those of a single SU(2) Ooguri model for the same boundary data (k,n). 
This is because critical points of the gravity model of the form X + = X~ just correspond to 
critical points of the SU(2) 15j symbol, and, moreover, the action of the gravity models ([28)) 
reduces to the SU(2) action (|15p on setting e = e'. The Hessians may differ, however. In a sense 
this explains all of the asymptotic terms except the genuine gravity terms (the genuine gravity 
terms being those with asymptotics given by 7 times the Regge action, for which X + and X~ are 
inequivalent). A similar conclusion holds for the Lorentzian model based on the representations 
of SL(2,C). Apart from the expected gravity terms given by 7 times the Lorentzian Regge 
action, all of the other asymptotic terms correspond to the SU(2) Ooguri model. 

This picture can be understood heuristically from the Lagrangian point of view. The gravity 
models are obtained by constraining a BF theory and so the Lagrangian should be a general- 
isation of the Reisenberger Lagrangian |Rei981 lBHN + 04] . which corresponds to our case with 
c 7 = 1. The Reisenberger Lagrangian is the SO(4) or SO(3, 1) BF Lagrangian with a constraint 
term added which enforces the equality of the area metric from the self-dual part with the area 
metric from the anti-self-dual part of the B field. One class of solutions is given by 

n IJ — l J J C K a C L 

B ~ 2 KL e A e ' 

which correspond to our 'geometric' boundary data and the gravity asymptotic terms. A second 
class of solutions is given by an SU(2) connection A and the matrix B IJ lying in the su(2) Lie 
algebra, for which the action is the SU(2) BF action with no constraints. This class of solutions 
corresponds to our non-geometric ones. 

Finally, some remarks about the full state sum models and their critical points. The question 
of whether the Regge equations of motion appear correctly from the state sum models is an 
important question. From the results developed here, it is possible to draw some preliminary 
conclusions. In the case of the Euclidean models, geometric boundary data for a 4-simplex is 
always embedded in a larger space of non-geometric boundary data. This can be seen from the 
fact that the vectors b a f, span an 18-dimensional vector space but the geometric ones form a 
submanifold of dimension only 13. Thus one can always perturb the h a i, away from the geometric 
confiurations, and hence one can perturb the boundary data given by the away from the 
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geometric configuations too. Thus at first sight it seems that the variations are more than the 
variations of the metric geometry in Regge calculus, and hence one will get more equations of 
motion. However the fact that the number of solutions drops immediately on leaving the space 
of geometric solutions means that the situation is a little more complicated; there is interference 
between the geometric sector and the SU(2) BF sector. If one could separate out the geometric 
solutions, then these would satisfy the Regge equations, but the SU(2) BF ones certainly do 
not. More work is required to determine the status of the Regge equations in this model. 

For the Lorentzian theory however, the situation is much simpler and clearer. The gravity 
solutions are well separated from the SU(2) BF solutions and all though the state sum model 
contains both in its asymptotics, perturbing one does not give the other. Thus the Lorentzian 
gravity configurations can only be perturbed to other Lorentzian gravity configurations. This 
means that the correct Regge-like gluing of 4-simplexes is preserved under perturbations. Thus 
perturbations can lead only to the Regge equations of motion. Thus the state sum model does 
indeed contain asymptotic terms which approximate general relativity. 
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